Low temperature properties of core-softened models: water vs. silica behavior 
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A core-softened model of a glass forming fluid is numerically studied in the limit of very low tem- 
peratures. The model shows two qualitatively different behaviors depending on the strength of the 
attraction between particles. For no or low attraction, the changes of density as a function of pres- 
sure are smooth, although hysteretic due to mechanical metastabilities. For larger attraction, sudden 
changes of density upon compressing and decompressing occur. This global mechanical instability 
is correlated to the existence of a thermodynamic first-order amorphous-amorphous transition. The 
two different behaviors obtained correspond qualitatively to the different phenomenology observed 
in silica and water. 
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I. INTRODUCTION 

The recently acknowledged possibility of the existence 
of single component systems which display coexistence 
between two different liquid phases has opened many in- 
teresting questions, and shed new light into thej-sUidy of 
the anomalous properties these systems displaycTTJ. The 
case of water is probably the most intensively studied, 
due to its ubiquity in nature. There is by now a general 
consensus that water displays a transition between two 
different amorphous states in the supercooled region of 
its phase diagramo. Experiments carried out in water at 
T ~ 130^ show an abrupt change of volume v as a func- 
tion of pressure P, which indicates the existence of the 
first order transitionQ. The v(P) curve is hysteretic, and 
the jump in v occurs at P ~ 0.3 GPa upon compressing, 
and at P ~ 0.05 GPa upon decompressing. The volume 
change at the transition is about 0.2 cm 3 /g. The two 
amorphous phases of water are smoothly related to two 
different liquid phases at higher temperatures, which at 
coexistence determine a first order liquid-liquid transition 
line ending in a critical point, located in the metastable 
region of the phase diagram. Most of the anomalies of 
water are usually interpreted as a consequence of the ex- 
istence of this liquid-liquid line and the liquid-liquid crit- 
ical point, but in general the existence of liquid-liquid 
coexistence is not a necessary condition for the existence 
of most of the anomalies!! 

Water is not an isolated case. There is a whole family 
of substances, usually referred to as tetrahedrally coor- 
dinated materials, that display many of the anomalies of 
wateitl. Within this family, another particularly inter- 
esting case is that of amorphous SiC>2 (silica). Many of 
the anomalous properties of water are also found in sil- 
ica. Some of them have been observed in experiments (as 
the density anomaly) and other only in numerical sim- 
ulations up to now (as maxima of isothermal compress- 
ibility and diffusivity as a function of pressures). This 
has led to think that possibly a first order amorphous- 
amorphous transition also occurs in silica. But evidence 



of this transition has been elusive. Experiments at am- 
bient temperature in silica show an irreversible increase 
of density when the system is compressed up to P > 20 



This behavior is 
However, for no 



GPa and successively decompressed 
reproduced in numerical simulations 
particular value of pressure there is a sudden change in 
density that could be interpreted as a direct evidence of 
a first order transition. 

It has been argued that the qualitatively different be- 
havior of water and silica is due to the temperature at 
which experiments are carried out .compared to the glass 
temperature T g of the materialMeEl Whereas T/T g is 
about 0.1 for the experiments in silica, it is close to 1 for 
water. Some people have raised the expectation that if 
compressing experiments were done in silica at tempera- 
tures near or above T g , they would reveal the first order 
transition, which is supposed to be hidden in the ambient 
temperature experiments due to lack of thermodynamic 
equilibrium. 

We will not study the actual behavior of water and 
silica. Instead, we will study a very simple model of a 
glass former at low temperatures, including the limiting 
T = case, in different regions of the parameters to gain 
qualitative insight into the real cases of water and silica. 
The model consists of particles interacting through a po- 
tential with a hard core plus a soft repulsive shoulder. 
In addition, an attractive contribution to the interaction 
can be included. This kind of model is not unrealistic for 
the study of the properties of tetrahedrally coordinated 
materials, and in fact it was shown to xenroduce many of 
the anomalies these materials possess&EJ. We will show 
that when there is no attraction between particles (or 
only a weak one) , implying in particular that there is no 
first order transition, the v(P) curve of our model in sim- 
ulations at T — shows an hysteretic behavior, which 
is associated to the existence of different mechanically 
stable configurations, and qualitatively agrees with the 
known phenomenology of silica. The inclusion of a suffi- 
ciently strong attraction may produce the appearance of 
a first order transition, which is clearly observed even in 
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simulations at T = 0, in the form of a global mechanical 
instability. The form of v(P) we obtain in this case is 
very similar to what is found in water. 

These findings will lead us to propose a different sce- 
nario to place together the properties of water and sil- 
ica. We suggest that the first order transition in water 
will be observable in the form of a mechanical instability 
even in experiments at T — > 0. On the other hand we 
stress the possibility that there is no first order transi- 
tion in silica, and that experiments performed in silica 
near T g will not reveal new qualitative ingredients. It is 
worth noticing that the existence of the thermodynamic 
anomalies of silica does not contradict this interpreta- 
tion, since a liquid-liquid critical point is not necessary 
to observe those anomalies, as it is knownJrom studies 
in models closely related to the present oneS. 

The paper is organized as follows. The model and de- 
tails on the simulation procedure are presented in Section 
II. The results are contained in Section III, and the rele- 
vance of them to silica and water is presented in Section 
IV. Section V contains the summary and some final com- 
ments. 



II. THE MODEL 

The model we study is defined as follows. We consider 
a bidisperse set of spherical particles, in order to avoid 
crystallization. Particle i is characterized by the value 
of a parameter n, which is taken from a bimodal distri- 
bution, i.e., Ti = r a , or fj = Tb, with equal probability. 
The interaction potential U between particles i and j de- 
pends on f = r/R where r is the real distance between 
particles, and R = + rj . The potential U is composed 
of a repulsive and an attractive part, U = U R + U A . The 
form we use for U R is 

U R = 00 for f < 1 

U R = e - f— + 1.2-1.8(f-l.l) 2 V for Kf< 
r \r — 1 / 

U R = for f > 1.9202. 

where r — r a + rv U R is plotted in Fig. [l| The form 
of this potential is a smooth version,pf. a potential that 
we have studied in detail previouslytrE-3. This smoother 
form is preferred here in order to avoid ambiguities in 
the calculations at T = 0, that appear in case the forces 
are not continuous. The particular analytical form we 
use is not really important, the only crucial feature of 
the potential is the existence of two preferred distances 
between particles (ro ~ 1.1 R and r\ ~ 1.9 R in Fig. 
|l|) depending on the value of the external pressure. The 
attractive part of the potential U is simply given by 

R 

U A = s —a(r - b) for f <b 

r 

U A = for r > b (2) 



with the two dimensionless parameters b and a fixing re- 
spectively the range and intensity of the attraction. Two 
examples of the potential with the attraction term (those 
to be used in the simulations) are also shown in Fig. [l]. 
In all the results to be presented (corresponding to two- 
dimensional systems), we use r a = 0.45f, and r\, — 0.55f. 
Temperature will be measured in units of eok^ 1 , pressure 
in units of eo/c^r -2 , and volume in units of r 2 . 

We simulate the system by standard Monte Carlo tech- 
niques. Particles are placed in a box with periodic bound- 
ary conditions. At each time step the position of a single 
particle is modified to a new position which is randomly 
chosen within a sphere of radius 0.01 f centered at the 
original position. This trial movement is accepted ac- 
cording to the Metropolis rule. The update is made se- 
quentially for all particles. Results using constant pres- 
sure simulations and others at constant volume will be 
shown. In the constant pressure scheme, the volume of 
the system is considered as an additional Monte Carlo 
variable, and homogeneous expansion and contraction of 
the coordinates of all particles (and also of the size of 
the simulation box) are tried. This permits the volume 
of the system to adjust to the given external pressure. 
In constant volume simulations, pressure is calculated as 
minus the energy change divided the volume change in a 
small (virtual) homogeneous rescaling of all coordinates 
of the particles. 

III. RESULTS 

We will show results for a two dimensional system. 
This is done to allow a better visualization of the con- 
figurations, however we emphasize that all the results 
we discuss were qualitatively re-obtained in three dimen- 
sional systems. We start by showing results in the case 
of no attraction (a = 0), at T = 0. The particles are 
randomly placed in space at the beginning of the simula- 
tion, and the system is rapidly compressed up to reaching 
a mechanically stable configuration at P ~ 0.5, v ~ 3.2. 
Then P is increased, or v is decreased (depending on the 
kmd of simulations) by small steps. Mechanical equilib- 
rium is obtained at each P-v value. In this way we reach 
P ~ 4, v ~ 1.5. Then we slowly expanded the system to 
the original values of P and v. 

We see in Fig. || the values of P and v during this 
process, both in simulations at constant P an at con- 
stant v for a system of 200 particles. The curves show 
a series of small though abrupt changes in v or P (de- 
pending on the kind of simulations) that correspond to 
mechanical instabilities in the system. Also the hystere- 
sis in v(P) we observe (which we found is repetitive upon 
compressing and decompressing, though only one cycle is 
shown in Fig. ||) originates in the existence of mechan- 
ical metastabilities in the system. An indication of this 
fact can be obtained from the following arguments. Let 
us define v = dh/dP, where h = e + Pv is the enthalpy 
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per particle of the system as obtained from the simula- 
tions. If the system was in thermodynamic equilibrium 
we should obtain v = v, since this is a thermodynamic 
identity at T = 0. We plot v (from the constant P sim- 
ulations) also in Fig |2[ We see that there are systematic 
differences between v and v, particularly in the region 
where v changes rapidly. The difference between v and v 
is due to mechanical instabilities upon compressing and 
decompressing that produce energy to be dissipated in 
the process. A very simple example is illuminating in 
this respect. For two particles interacting with the po- 
tential of Fig. [I] (continuous line) at T — 0, the evolution 
of distance d (that replaces volume in this case) and en- 
thalpy h as a function of the compressing force F (that 
takes the role of pressure) is shown in Fig. [| The me- 
chanical hysteresis in d(F) is precisely the same effect 
we observe in the simulations, the only difference is that 
in the simulations the disordered nature of the system 
produces a smoothing of v(P). In the two particle prob- 
lem d = dh/dF equates d except by the existence of 
two delta peaks at the points where there is a jump in 
h. This is the energy being dissipated. In the numeri- 
cal simulations there is an averaging over many different 
atomic environments, and the delta peaks are smoothed, 
but still clearly visible as regions where v < v during 
compressing, or v > v during decompressing in Fig. |[ 

The mechanical instabilities in the system are of a local 
nature, in the sense that they produce non-correlated re- 
arrangements of particles as pressure or volume changes. 
Then we expect the small reentrances in the constant v 
simulations that signal the existence of this instabilities 
to become weaker in larger samples, since they are aver- 
aged over the whole system. In fact, as an indication of 
this, we show in Fig. [| results for a system of 1000 parti- 
cles, as compared to the case of 200 particles, and we see 
that 'fluctuations' are considerably smaller. The global 
amount of hysteresis however remains quite the same. 
Fig. H shows snapshots of the 1000 particles system at 
the points indicated in Fig. Qand Fig. ^ shows the cor- 
responding radial distribution function. We see that as 
the pressure increases and volume decreases, more and 
more particles move from a typical distance n ~ 1.9R 
between neighbors to the shorter distance tq ~ l.liZ. 
The collapse of neighbor particles from r\ to tq is not 
a collective effect, it occurs in a non-correlated way in 
different positions of the sample, namely no abrupt tran- 
sition exists. 

From the two particle problem we see an interesting 
characteristic of the h(F) curve. The values of h in 
the compressing and decompressing branches cross each 
other (at Ff/eo ~ 2.3). This crossing is also observed in 
the complete simulations. In fact, in Fig. [?] we see values 
of h obtained during compression and decompression, in 
the simulations at constant P (all results are for a system 
of 1000 particles, from now on). We see clearly a cross- 
ing of the two branches at P ~ 1.8. This crossing has no 
profound physical meaning, and in particular it docs not 



indicate the existence of a first order transition between 
compressing and decompressing branches, since these are 
not the thermodynamic values of the free energy, and the 
states accessed during compression and decompression 
do not exhaust the whole phase space of the system. To 
show this in more detail, we made simulations in which 
the system was annealed at fixed volume, trying to get as 
closer as possible to the real ground state of the system at 
each P. The final values of v that we get as T — > at dif- 
ferent P are shown in Fig. ||. Although we cannot guar- 
antee that the true ground state was obtained, the values 
of h(P) in this case are systematically lower than both 
the compression and decompression ones, as we can see 
in Fig 0. For the values of v{P) obtained in the annealing 
simulations, the thermodynamic relation v = dh/dP is 
well satisfied within the numerical errors (Fig. |^). The 
values of v(P) obtained in the annealing process are in- 
side the loop of compression-decompression, as it could 
be expected for the equilibrium values, since the hystere- 
sis loop represents the maximum amount of mechanical 
metastability that the system is able to sustain. In addi- 
tion, v(P) is continuous, and this is consistent with the 
lack of a thermodynamic first order transition at T = 
in our model. 

To obtain a thermodynamic first order transition, and 
to study the way in which it reflects in the compres- 
sion decompression results at T = 0, a certain amount 
of attraction between particles must be included. The 
existence of attraction in the system may have profound 
effects in its phase behavior. We will first discuss ana- 
lytically the case in which an attraction energy 7(7') < 
of infinite range is added to the interparticle energy U. 
In this case, the energy per particle e gets an additional 
contribution Se of the form Se = J v S(r)"f(r)dr, S(r) be- 
ing the radial distribution function. The assumed limit 
in which the range of j(r) goes to infinity, corresponds 
to the case in which this integral is governed only by the 
limiting value of S(r) as r — > 00, i.e., it will depend only 
on the density of the system. In this case, S(r) itself is 
not affected by the attraction term. Then Se takes the 
form Se = —j/v with some constant 7 > 0. This term is 
directly added to the free energy G of the system. When 
computing the equation of state from dG/dv — 0, the 
only difference with the case without attraction is that 
P in the state equation is replaced by P + 7/w 2 . Then 
from the results of the simulations without attraction we 
can immediately get those for a system with an attraction 
that has infinite range, just rescaling (self-consistently) 
the pressure axis through P — » P + 7/w 2 . The result 
of this procedure is shown in Fig. Q We see that due 
to the volume-dependent rescaling of P, and to the form 
of the v(P) curve with no attraction, a re-entrance in 
v(P) may appear if the attraction is higher than some 
minimum. The reentrance in the thermodynamic v(P) 
curve indicates the existence of a first order transition be- 
tween two amorphous phases of different densities. This 
reentrance appears also for the compressing and decom- 
pressing branches, basically at the same value of 7. This 
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means that in this case the mechanical instability at zero 
temperature exist if and only if the system has a thermo- 
dynamic first order transition at T = 0. Then the global 
mechanical instability upon compression and decompres- 
sion as pressure passes through the spinodal points at 
which dv/dp becomes infinity (points C and C in Fig. 
^) is an indirect observation of the first order transition. 
We emphasize that in the present case the instability 
is global, in the sense that once it occurs, it involves a 
finite fraction of the whole system, contrary to instabili- 
ties in the case without attraction, which are associated 
to individual particles. We note that for a truly infinite 
range attraction (more precisely, if the range of the at- 
traction is much larger than the system size) the loop 
in v(P) when there is a first order transition is physical, 
namely, it is observable in constant volume simulations, 
and no Maxwell construction can be invoked to flatten 
it out. In constant P simulations instead, we would get 
an abrupt volume change when we go through the pres- 
sure corresponding to the spinodal points, and this is 
what we are referring to as a global mechanical instabil- 
ity. From Fig. || we also see that the position of the first 
order loop moves towards lower pressures as the strength 
of the attraction is increased, and then it can be com- 
pletely moved into the metastable P < region when 
the attraction is strong enough. Notice also that in cases 
in which the compressibility anomaly of the purely re- 
pulsive model (the rapid change in volume around P ~ 2 
in our case) is weaker, it may happen that no first order 
transition appears at all, for any value of the attraction. 

The previous analysis tells us that in the case of a 
long range attraction, the existence of a first order ther- 
modynamic transition and mechanical instabilities in the 
compression-decompression path at T = are closely re- 
lated, each of them implying the existence of the other. 
We want to analyze now to what extent this scenario can 
be extended to the case in which the attraction is short 
ranged. Then we conducted compressing-decompressing 
simulations using the finite range form (||) of the at- 
traction. We first use a rather large value b — 4 for 
the attraction range, expecting to reobtain basically the 
mean field phenomenology. For weak attraction [c < 0.1 
in (||)], we get v(P) curves which are slightly modified 
with respect to the case of no attraction, but no qualita- 
tively new results appear. But if the attraction is strong 
enough, we obtain signs of global mechanical instabilities 
and a first order transition. The results of compression- 
decompression simulations are illustrated in Fig. ^ for 
c = 0.15. We see that constant v simulations get a re- 
gion (for 1.8 < v < 2.4 upon compression) with small 
reentrances in the calculated values of P. In turn, in 
constant P simulations we see an abrupt collapse of vol- 
ume that jumps in a finite amount. This jump signals 
the occurrence of a global mechanical instability. 

The snapshots of the system when going through the 
coexistence region (Fig. |ll|) show indeed that bubbles of 
the two amorphous phases coexist. Note the difference 
between the configurations in Fig. O, and those with no 



attraction during compression and decompression [Fig. 
H(B) and (D)]. Here the particles that collapse to the new 
phase tend to form well defined clusters in the sample, 
whereas in the other case the system remains uniform. 
This collective clustering is the driving force of the me- 
chanical instability. 

In this case in which the attraction range is rather large 
we expect the mean field arguments of the previous sec- 
tion to apply, and then that the system has (in addition 
to the global mechanical instabilities at zero tempera- 
ture) a thermodynamic phase transition at sufficiently 
low temperatures. To check directly for this transition 
is of course not an easy task in general, but for the cur- 
rent parameters it turns out that the transition is still 
observable in equilibrium simulations at finite tempera- 
ture, and then we can be sure that it will persist down 
to zero temperature. In fact, to show only one indicator, 
the v(P) function in constant v simulations at T — 0.09 is 
also plotted in Fig. |l^. We see that hysteresis upon com- 
pressing and decompressing has completely disappeared, 
as it should be since we are in thermodynamic equilib- 
rium. The v(P) curve has an abrupt change at P ~ 0.75, 
signaling the existence of a two phase coexistence region, 
and then a first order transition. 

Next we simulated a case in which the attraction is 
much shorter ranged, specifically, we used the attraction 
form (^) with a = 2.2 and c = 1.0 (the interaction po- 
tential in this case can be seen in Fig. [j]). For these 
parameters we are still able to detect the existence of 
the first order thermodynamic transition in simulations 
at finite temperature. In Fig. [l^ we plot the v(P) curve 
at T — 0.09, in which hysteresis has disappeared, and 
signs of a first order transition at P ~ 0.55 are appar- 
ent. The compression-decompression curve at T = is 
shown also in Fig. [l2| We see that in this case in which 
the attraction is much shorter ranged, the signs of global 
mechanical instabilities are still clearly observable. Note 
also that in this case, the decompression to P = (in the 
constant P simulations) was such that the system was 
not able to regain its initial density, and remains in a 
densified structure. 

In simulations performed with attraction of even 
shorter range, or when this is not strong enough, we did 
not observe global mechanical instabilities, but neither 
a thermodynamic first order transition, at least in the 
temperature ranges where equilibrium calculations can 
be done. 

Let us briefly summarize the numerical results we 
have presented, before going to their relevance to un- 
derstand the phenomenology of water and silica. For the 
case of no attraction between particles the v(P) curve 
at zero temperature shows hysteresis upon compression- 
decompression, which is originated in local mechanical 
metastabilities. The thermodynamic v(P) in this case is 
smooth and globally stable (dv/dP < 0) When an at- 
traction of sufficiently long range is included, and in the 
case this is stronger than some minimum value we get 
both a thermodynamic phase transition and global me- 
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chanical instabilities at T = 0. For attractions of shorter 
range, and in cases we can guarantee the existence of a 
first order transition, we have observed that the global 
mechanical instability still exists. 

IV. UNIFYING THE SILICA AND WATER 
PHENOMENOLOGY 

Silica and water display many thermodynamic anoma- 
lies. They include the well known density anomaly (a 
temperature at which density is maximum), compress- 
ibility and specific heat anomalies, and diffusivity anoma- 
lies. These anomalies do not require the existence of a 
first order liquid-liquid or amorphous-amorphous transi- 
tion to exist. In fact, basically all of them are found in 
the present model with no attraction, i.e., wiien we can 
be sure that there is no first order transitionB. 

Experimental and numerical evidence suggests that 
water possess a first order liquid-liquid equilibrium at 
low temperatures. Experimental evidence of a first or- 
der transition in silica is lacking. The present model 
shows that hysteresis in v(P) curves is by no means a 
strong evidence of a first order transition if a global me- 
chanical instability is not observed. And in fact, exper- 
iments show that silica is globally stable upon compres- 
sion and decompression, although displaying hysteretic 
behavior. Some numerical evidence of a first order trait, 
sition in silica has been presented. In one case howeverE3 
the evidence was just the crossing of the free energies 
obtained during compression and decompression simula- 
tions at T — of model silica. This was erroneously 
attributed to an underlying first order transition. As we 
showed in our model this crossing occurs whether there 
is a first order transition or not, and it is due to micro- 
scopic metastabilities. More seppus evidence come from 
the simulations by Voivod et aZ.EJ. They use two different 
numerical models of silica and extrapolate high temper- 
ature results to zero temperature, finding evidence of a 
first order transition. Still, the results are preliminary, 
and based on extrapolations that call for more detailed 
study. 

We have shown that a very simple and transparent 
model has or has not a first order transition depending on 
the strength of the attraction that is included. We have 
shown simulations in which a mechanical instability at 
T = and a thermodynamic first order transition exist. 
The mechanical instability is qualitatively similar to that 
found in water in experiments at T ~ T g . This leads us to 
expect that water must display the same instability even 
in experiments at much lower temperatures. To confirm 
this expectation it would be interesting if some of the 
models that are used to simulate water, and that have 
shown liquid-liquid coexistence, were tried in the limit of 
zero temperature to look for the mechanical instability, 
or even if compression decompression experiments were 
done in water at much lower temperatures. 



In the case of a long range attraction between parti- 
cles, we have shown that a thermodynamical first order 
transition and global mechanical instabilities are closely 
correlated. It is tempting to expect that for short range 
attraction this correlation remains, but this is not guar- 
anteed. In cases in which attraction is short ranged, and 
when the first order thermodynamical transition is ob- 
servable in equilibrium simulations, we have still observed 
global mechanical instabilities. We have not been able 
to identify a set of parameters for our model where a 
thermodynamical transition occurs in the lack of global 
mechanical instabilities. But this might be due to our 
incapacity of detecting a thermodynamic first order tran- 
sition at very low temperatures when this transition does 
not show up at temperatures at which equilibrium simu- 
lations are possible. 

On the basis of the present results we consider that 
the two possible scenarios for silica are the following. It 
may happen that silica does not have a first order tran- 
sition. Our results for zero or low intensity attraction 
correspond in fact to a case in which a first order tran- 
sition does not exist, and the phenomenology we obtain 
is closely related to that of silicaE-J Still, the other open 
possibility is in fact that silica has a first order tran- 
sition but this is not reflected in the zero temperature 
compression-decompression experiments. In this respect 
again, it would be nice if the numerical models used for 
silica that seem to posses a first order transition (and 
in case in fact this is confirmed) were used at T = 
to search for a global mechanical instability. If this is 
obtained, then those models will be shown to be no re- 
liable to describe silica in this limit, since we know real 
silica does not have this instability. If the model does not 
show a mechanical instability instead, this would indicate 
that the correlation between mechanical instabilities and 
a first order thermodynamical transition is not universal. 

V. SUMMARY AND CONCLUSIONS 

In this work we have studied a model of a glass form- 
ing fluid, consisting of spherical particles with a hard 
core, a repulsive shoulder and some amount of attrac- 
tion. We investigated to what extent the existence of 
global mechanical instabilities at T = is correlated to 
the existence of a thermodynamical first order transition 
between two amorphous phases. We have obtained that 
the model without, or with weak attraction between par- 
ticles does not display neither a first order transition nor 
global mechanical instabilities. This behavior coincides 
with the known phenomenology of silica. In the presence 
of a strong attraction between particles, both a first order 
amorphous-amorphous transition and global mechanical 
instabilities at T = were obtained. This scenario cor- 
responds to the behavior of water. Our work suggests 
that the difference between the phenomenology of silica 
and water may be related to the lack of an amorphous- 
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amorphous transition in silica, in opposition to the exis- 
tence of this transition in water. This is not incompatible 
with the existence of coincident thermodynamic anoma- 
lies in both cases. 
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FIG. 1. The interparticle potential U(r) without attraction 
(continuous line) and for two different choices of the attractive 
part. 
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FIG. 2. v(P) function for the system without attraction, at 
T = 0. Results of simulations at constant P and constant v 
are shown. Also shown are the values of V = dh/dP, from the 
constant P simulations. Differences between v and v reveal 
lack of thermodynamic equilibrium. 
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FIG. 3. Distance d (a) and enthalpy h = U + Fd (b) for two 
particles interacting with the potential of Fig. [l] (continuous 
line) as a function of the compressing force F. Compressing 
and decompressing routes are indicates. 
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FIG. 4. v(P) curve as in Fig. ^| for systems with 200 and 
1000 particles. Fluctuations tend to average out in the larger 
sample, but the overall hysteresis remains the same as in the 
case with fewer particles. 
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FIG. 5. Snapshots of the system at the points indicated in 
Fig. 0. Dot size represents the hard core of the particles. 
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FIG. 6. Radial distribution functions of the configurations 
shown in Fig. ^. The weight transfer between ro ~ 1.17? and 
n ~ 1.97? is clearly visible (the triplet structure of the peaks 
originates in the bidispersity of the system). 
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FIG. 7. Evolution of enthalpy h as a function of pres- 
sure during compression and decompression, and from in- 
dividual simulations annealing the system at fixed pressure. 
Compressing and decompressing branches cross each other at 
P ~ 1.8, but the thermodynamic values are always lower. 
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FIG. 8. v(p) curve obtained as the limit T — > of individ- 
ual simulations at constant volume (large full symbols), and 
calculated v = dh/dP from the h values obtained in the same 
simulations (open symbols). Within the numerical errors we 
get v = v, is it should be in thermodynamic equilibrium. The 
compressing-decompressing loop of Fig. |i| is also shown to 
allow comparison. 




FIG. 9. The effect of an infinite range attraction of differ- 
ent intensities on the v(P) curves of Fig. |^. A reentrance, 
indicative of a first order transition is observed for the two 
largest values of 7. For 7 ~ 9eor 2 the first order transition is 
moved into the metastable P < region. Large symbols are 
the (almost) thermodynamic values, and small symbols are 
the compressing and decompression results. The position of 
the spinodal points C and C' upon compression and decom- 
pression is also shown in the case 7 = 7eof 2 . 
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FIG. 10. Full simulations with an attraction term of the 
form (Q) with b — 4, and a — 0.15. Constant pressure and 
constant volume simulations are shown. Note the abrupt v 
change in the constant P simulations, and the weak reen- 
trances of the constant v simulations. These are mechan- 
ical instabilities that indicate the existence of a thermody- 
namic first order amorphous-amorphous transition in this 
case, which is made evident in the equilibrium v(P) curve 
at T = 0.09. 




FIG. 11. Snapshots of the system at the points indicated 
in Fig. [lo| Note the coexistence of rather large clusters of 
two different phases with different densities, as for instance in 
the encircled regions (compare with Fig. 
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